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The transport properties of organic light-emitting diodes in which the emissive layer is composed 
of conjugated polymers in the liquid-crystalline phase have been investigated. We have performed 
simulations of the current transient response to an illumination pulse via the Monte Carlo approach, 
and from the transit times we have extracted the mobility of the charge carriers as a function of 
both the electric field and the temperature. The transport properties of such films are different from 
their disordered counterparts, with charge carrier mobilities exhibiting only a weak dependence on 
both the electric field and temperature. We show that for spatially ordered polymer films, this 
weak dependence arises for thermal energy being comparable to the energetic disorder, due to the 
combined effect of the electrostatic and thermal energies. The inclusion of spatial disorder, on the 
other hand, does not alter the qualitative behaviour of the mobility, but results in decreasing its 
absolute value. 
PACS: 72.10.-d, 72.80.Le 



I. INTRODUCTION 

Conjugated polymers, which derive their semiconduct- 
ing properties from delocalized ir bonding along the poly- 
mer chair£l, enable cheap, lightweight, portable and large 
area displays, with jj^omising applications to solar cells 
and photodetectorsBlm Since the ficst observation of 
polymer light-emitting diodes (LEDs)B, detailed inves- 
tigation of the aspects of chemistry, physics and [Engi- 
neering of materials has resulted in rapid progressu. In 
order to commercialize polymer devices, high efficien- 
cies, brightness, and carrier lifetimes are required!! It is, 
therefore, essential to fully understand the fundamental 
physics of electrical transport through conjugated poly- 
mers, for a recent review seeo. 

Previous theoretical studies have dealt with strongly 
disordered organic materials, in^zhidp charge transport 
is mainly attributed to hoppingoEJO. The charge car- 
rier mobilities of such systems are low and exhibit a 
strong temperature and electric field dependence. How- 
ever, large carrier mobilities arc generally highly de- 
sirable-, [Most recently, measurements in light-emitting 
diodealj : E£lEJ have demonstrated enhanced carrier mo- 
bilities, varying only weakly with the electric field, for 
devices characterized by high degree of order in the poly- 
meric material. This weak dependence, which contradicts 
the theoretical predictions, has been attributed to the pu- 
rity of the polymer films, but it is not yet understood. 

This realization has motivated us to investigate 
the transport properties of liquid-crystalline conjugated 
polymer films in which the chains are nematically aligned 
perpendicular to the direction of transport. In a pre- 
vious worklij we looked into the character of transport 
through such polymer films. By employing the Monte 



Carlo technique, we showed that it is possible to obtain 
non-dispersive transport in such systems. This observa- 
tion is in agreement with time-of-flight exper imen ts con- 
ducted on liquid-crystalline polyfluorene filmsE3lI3, which 
demonstrate non-dispersive hole transport with enhanced 
charge carrier mobilities compared to previously exam- 
ined conjugated polymers. Moreover, we established 
the conditions under which such transport is retained. 
However, our initial model was a simplistic one, with 
the film morphology reduced to occupied sites of a two- 
dimensional lattice, the variation of the hopping prob- 
ability with the field being rather crude, and there was 
no temperature dependence. In order to accomplish more 
detailed studies, we have first extended the model to fully 
describe the geometry of the polymer film, so that various 
film morphologies, similar to those appearing in realistic 
systems, can be considered. Secondly, we have explicitly 
included all parameters of interest, such as the electric 
field, temperature, and disorder (both spatial and ener- 
getic). 

The aim of the present work is to probe in detail 
the field and temperature dependence of the mobility in 
liquid-crystalline conjugated polymer films. In particu- 
lar, we discuss the effect of the electric field on the inter- 
chain mobility of such polymer films, and we attempt to 
explain the weak dependence on the field. The interplay 
between electric field and temperature on the transport 
characteristics is also examined. We begin our investi- 
gation with spatially ordered films in which all polymer 
chains are perfectly aligned perpendicular to the direc- 
tion of the field. The additional effect of spatial disorder 
in the film configurations on the transport properties of 
such films, is also considered. To account for the chemi- 
cal regularity and the extended backbone conjugation of 



liquid-crystalline chains we have included only a small 
amount of energetic disorder. 



II. THE MODEL 

In order to investigate the transport properties of 
charge carriers moving within a liquid-crystalline poly- 
mer film, we have performed numerical simulations of the 
time-of-fiight technique. The system under study con- 
tains a polymer film sandwiched between two electrodes, 
and charge carriers are generated on one side by illumi- 
nation of the electrode with an intense pulse of light of 
short duration. The photogenerated carriers move within 
the bulk of the film under the effect of an external bias, 
whose sign determines the type of the charge carriers that 
will get transported. 

The film is composed of conjugated polymer chains of 
length L = lOOnm, which are nematically aligned per- 
pendicular to the direction of the electric field, chosen as 
the x direction. For the construction of the film periodic 
boundary conditions have been applied along the other 
two directions. Based on the extended backbone conjuga- 
tion of liquid-crystalline polymers, such as polyfluorene, 
that makes them stiff, and on bond vibrations being of 
very high frequency and low amplitude, we have assumed 
that the polymer chains can be described as rigid rods. 
The thickness of the film in the direction of transport has 
been taken equal to d = 1/zm, since relatively thick films 
are required for such measurements. 

Hopping motion of the charge carriers under the in- 
fluence of the electric field is assumed, which in general 
can be either intra- or inter-chain. This is described by 
the Monte Carlo technique, in which the carriers perform 
random walks within the bulk of the film until they reach 
the collecting electrode. The probability of hopping be- 
tween two sites i and j is equal to 



Vij — Pij j Pi j 



(1) 



where is the unnormalized hopping probability, which 
is taken to be of the form 



Pij = 7exp(- 



eE- 



2K B T 



(2) 



E is the electric field, ry is the relative position vec- 
tor, Kb is the Boltzmann factor, T is the temperature, 
and 7 denotes the electronic wavef unction overlap. Here, 
we have considered only nearest-neighbouring hopping 
within a cutoff distance equal to lOA. is the on-site 
energy of site i and is taken from a Gaussian distribution 



Pie) 



\/2n<r 2 



exp{-e 2 /2a 2 ) 



(3) 



whose width a determines the degree of energetic disor- 
der present in the film. Notice that in Equation |^ there 
is no other activation energy except the difference in the 
electronic site energies a carrier has to overcome during 
each jump. Equation ^| satisfies the principle of detailed 
balanceM 

Intra-chain charge motion is, in general_a very rapid 
process compared with the inter-chain oneEl Moreover, 
it was recently shown that in liquid-crystalline poly- 
mer filmsiater-chain hopping accounts for non-dispersive 
transport!!!! Taking these into account, we have mainly 
focused on inter-chain transport, so that in what follows, 
each Monte Carlo step corresponds to the time required 
for an inter-chain jump. 

The drift of the photogenerated carriers under the 
external bias. results in a time-dependent current I{t), 
which readaiH 



m = sii 



dxp{x)(^- 



1) 



(4) 



d is the film thickness, and p(x) is the charge density in 
the direction of transport, integrated over the y and z di- 
rections. From the current transient the transit time tr, 
which is the time the carriers need to reach the discharg- 
ing electrode, is obtained. We then use this extracted 
transit time to get the mobility \x via the relation 



(5) 



Et T 



III. RESULTS AND DISCUSSION 

We present results for the charge carrier mobilities of 
liquid-crystalline conjugated polymer films, obtained by 
using the Monte Carlo model described in the previous 
section. Unless stated otherwise, all inter-chain distances 
are equal to 10A. The output of such calculations is the 
current transient, Eq. ^, which arises from the applica- 
tion of the external bias. In Fig. [l] we present a cur- 
rent transient from our numerical simulations for elec- 
tric field E = 3 x 10 5 V/cm at room temperature, and 
for disorder a — 2, where a is defined as a = ct/KbTq 
and To = 300K. The calculated current exhibits the typ- 
ical behaviour of a time-of-flight signal, with a distinc- 
tive plateau followedJpy a decaying tail, thus shows non- 
dispersive transported. From the current plot we obtain 
the transit time i-r, which is required for the calculation 
of the mobility, by using the current integration mode. 
The latter gives the transit time as the point at which the 
current has fallen to half its value in the plateau region. 
In Fig. [j] the transit time is indicated by the arrow. 

In most experimental time-of-flight signals, the current 
transients initially exhibit a spike, which rapidly falls into 
the plateau. However, in our results such a pronounced 
peak is missingjij. We attribute the lack of the initial 



spike to the low degree of disorder present in our sys- 
tem, and in particular near the injecting electrode. A 
demonstration of the effect of purity, in the initial peak 
is presented in Fig. ||. Curve (a) shows the current tran- 
sient for a film with no energetic disorder at all (Se = 0) , 
whereas curve (b) corresponds to the case of a = 2. 
Curve (a) has been shifted along the y axis for clarity. Af- 
ter the illumination of the electrode, the simulated pho- 
togenerated charge packet begins its motion within the 
film. If no disorder is present, the packet is adequately 
described by a delta function moving within the bulk of 
the film with constant mean velocity. If some disorder is 
introduced into the system, the propagating packet will 
become more extended (Gaussian), and the carriers will 
no longer move together. The resulting change in the 
charge profile produces a spike in the current transient, 
as shown in curve (b), even though the disorder that we 
have considered is too small to give a pronounced peak. 
Additionally, spatial irregularities neat-the injecting elec- 
trode would lead to much larger peakstj, similar to those 
observed in experiments. 

By identifying the transit times from the current tran- 
sients for various values of the external electric field and 
for different temperatures, we have calculated the mobil- 
ity by using Equation 0. Figure || shows one set of such 
results for a = 1. In this graph the logarithm of the 
mobility, ln(/Lt), has been plotted as a function of \J~E~ for 
temperatures in the range of 100K to 350K. In this Fig- 
ure we see that for low temperatures ln/u shows a more 
pronounced dependence on y/~E, with a slight increase at 
low fields, which switches to a decrease at larger values 
of the bias. As the temperature increases, and, thus, the 
thermal energy becomes comparable to the energetic dis- 
order, this variation disappears and the mobility changes 
only weakly with the field. The reduction of the mobil- 
ity at large electric fields and at all temperatures is the 
outcome of the saturation of the transit time, since the 
electrostatic energy is sufficient for the charge carriers to 
move fast enough and exit the film at the minimum num- 
ber of time steps required. However, we should point out 
that we are not aware of any experimental results which 
show such a decreasing mobility. This indicates that in 
reality there are not liquid-crystalline polymer films of 
such energetic purity, even if they are characterized by 
perfect alignment of the chains. 

In Figure Hthe logarithm of the mobility is shown as a 
function of y/E and for the same range of temperatures 
as previously, but for disorder a — 2. In contrast to 
Fig. |3[ In /i now exhibits a positive gradient at all fields. 
For small values of T the thermal energy is not adequate 
to enable the carriers to surmount the energy barriers, 
and the electrostatic energy becomes important, ln/z in- 
creases with the field as the electrostatic energy permits 
the carriers to follow shorter paths within the film. Nev- 
ertheless, for larger temperatures, the energetic disorder 
is rather small, and the thermal and electrostatic energies 
have similar contributions. Thus, the strong field depen- 
dence seen in more disordered systems is smoothed out 



and the carriers move with a constant velocity. A direct 
comparison with experimental datalij shows that our cal- 
culations for the mobility at room temperature and for 
a = 2 lead to the same qualitative behaviour. 

Figure || shows the behaviour of the mobility as a func- 
tion of temperature. The circles correspond to a = 1 and 
the squares to a — 2. In the former case, In/x decreases 
with T for both small (filled circles) and large (open 
circles) electric fields. Energetic disorder is relatively 
small so that the thermal energy dominates, making the 
charge carriers to remain longer within the film, following 
larger paths. Hence, the mobility decreases with temper- 
ature. This is analogous to carriers being scattered by 
phonons in conventional semiconductors. When a in- 
creases, on the other hand, we see that at a low elec- 
tric field the mobility increases with temperature (filled 
squares), whereas for a larger field ln^i decreases with T 
(open squares). For small fields, thermal activation is the 
leading mechanism for the carriers to overcome the en- 
ergy barriers. Nevertheless, the energetic disorder is still 
weak and at large fields longer paths due to the excess 
thermal energy prevail, as before. This is an important 
result because it shows how much the mobility is affected 
by the interplay between electrostatic and thermal en- 
ergy in the case of small energetic disorder. A similar 
behaviour of the mobility as a function of temperature 
has also been observed in other types of polymers within 
a Master equation approachEJ, although in that case the 
crossover appears to occur at the same field for all tem- 
peratures considered. Moreover, at these small amounts 
of energetic disorder changes in In /z with temperature 
are also weak. 

Nevertheless, this crossover from increasing to decreas- 
ing mobility with temperature is not observed if we de- 
part from the equidistant chains arrangement. In Fig. ^ 
we show In /i as a function of \f~E~ and for various tem- 
peratures, for a film in which the liquid-crystalline chains 
are still all nematically aligned, but at inter-chain dis- 
tances that vary randomly between 7A and 13A (irreg- 
ularly spaced chains). The disorder strength is a = 2, 
as in Figure [|. A line has been fitted to the simulation 
data points to guide the eye. Hopping is only allowed 
between nearest-neighbours within the same cutoff dis- 
tance as before (10A), discussed in the previous section. 
When the chains are irregularly spaced within the film 
(Fig. |^) the mobility is enhanced with temperature for 
all values of the electric field, as-was shown in recent mea- 
surements of polyfluorene filmsEll. This behaviour resem- 
bles that of strongly disordered polymer filmsO. However, 
in our case the increase of the mobility with temperature 
arises from spatial irregularities in the arrangement of the 
chains within the film. The distribution in the relative 
position vector leads to variations in the electrostatic en- 
ergy (Eq. ^|), so that thermal activation dominates, and 
the field dependence of ln/i remains weak for all temper- 
atures. 

Comparison between Figures |i| and ^ also shows that, 
even though in both regularly and irregularly spaced 



polymer films ln/i varies only slightly with the electric 
field, the introduction of some spatial disorder leads to 
the reduction of ln/i. In the latter case, the charge carri- 
ers require more time to reach the discharging electrode 
as some nearest-neighbouring chains appear in distances 
larger than 10A. Such jumps are prohibited, and the car- 
riers must find alternative paths that might be longer. 
This would lead to greater transit times and, therefore, 
to smaller mobilities, as seen from Eq. |5| 

When the spatial disorder is retained and the energetic 
disorder is increased, ln/i decreases even more. More- 
over, its qualitative behaviour approaches the disordered 
limit, with a stronger dependence of ln/i on the electric 
field. This is shown in Figure [?], in which ln/i is plot- 
ted as a function of \/~E and at room temperature for 
an irregularly spaced film, for two different strengths of 
energetic disorder. The diamonds correspond to a = 2, 
and the circles to a = 3. The lines are drawn as a guide 
to the eye. Not only are the charge carriers obliged to 
avoid large inter-chain distances, but also they have to 
surmount larger energy barriers. In this case, the charge 
carriers rely mostly on the electrostatic energy as a source 
for the energy required to overcome the harriers, similarly 
to previously examined disordered filmsa. 

From the above discussion it is understood that the 
weak field and temperature dependence of ln/i in liquid- 
crystalline conjugated polymer films originates from the 
low energetic disorder present in the film. The addi- 
tion of spatial disorder, by irregularly spacing the liquid- 
crystalline chains, does not affect this dependence, but 
only results in the reduction of the absolute value of the 
mobility. Nonetheless, it is possible that another form 
of spatial disorder might be present in such films, and 
this is the existence of chains which deviate from perfect 
nematic alignment. In other words, it is possible that 
during the formation of such polymer films, some chains 
might appear at an angle 9 with respect to the axis of 
alignment. We have examined this case by assigning in 
our simulations a random value of this angle 9 to each 
chain. The chains were also positioned in random inter- 
chain distances, as in Figures ^ and [?]. 

Figure || shows the outcome of such calculations, in 
which 9 varies randomly between 0° and 5°. The loga- 
rithm of the mobility is plotted as a function of y~E, for 
various temperatures, as before, with energetic disorder 
a = 2. The lines are drawn as a guide to the eye. This 
Figure shows that ln/i is only weakly dependent on the 
electric field, and increases with temperature, in a simi- 
lar fashion to the case of perfect alignment. Hence, the 
inclusion of an extra form of spatial disorder does not 
seem to influence the qualitative behaviour of the mo- 
bility, and stresses the importance of low energetic dis- 
order as the factor that determines the degree of depen- 
dence on the field and the temperature. Unfortunately, 
we cannot make any comparisons as to whether the mo- 
bility increases or decreases with respect to the perfect 
alignment case, since our simulations for the randomly 
oriented case have been performed for films with thick- 



ness d — 0.5/im. Although such film thickness might 
be considered small for real time-of- flight experiments, 
we believe that our calculations are able to represent the 
qualitative behaviour of the mobility for films which show 
deviations from perfect alignment. 



IV. SUMMARY AND CONCLUSIONS 

Liquid-crystalline polymer films are generally expected 
to be characterized by spatial purity and chemical reg- 
ularity. In this paper, we have presented calculations of 
the inter-chain mobility of charge carriers through such 
polymer films, by employing the Monte Carlo technique. 
In particular, we have investigated the effect of low ener- 
getic disorder in the electric field and temperature depen- 
dence of the mobility, in order to account for the weak 
dependence observed in experiments. Initially, films with 
no spatial disorder were examined. Our numerical sim- 
ulations have shown that for relatively clean systems, 
a = 1, the mobility decreases with the field, as a re- 
sult of the saturation of the transit time by the elec- 
trostatic energy. Nevertheless, this behaviour does not 
depict a realistic situation and a higher amount of dis- 
order, a — 2, was required lo achieve qualitative agree- 
ment with experimentsEjO'El. When the thermal energy 
is not enough for the charge carriers to overcome the en- 
ergy barriers, the electrostatic energy becomes important 
and ln/i increases with the field. The weak field depen- 
dence of the mobility arises from equivalent contributions 
from the thermal and the electrostatic energy when the 
thermal energy is comparable to the energetic disorder. 

The interplay between energetic disorder and thermal 
energy also appears in the temperature dependence of the 
mobility at a given electric field. For low energetic disor- 
der, a = 1, the mobility decreases with temperature at 
all fields, as is commonly seen in crystalline conductors. 
For larger energetic disorder, a — 2, and for constant 
inter-chain distances, the mobility increases with tem- 
perature for small fields, but a crossover to the opposite 
behaviour appears at larger fields. This result is another 
demonstration of the combined effect of the electric field 
and the temperature on the mobility of charge carriers 
through spatially ordered liquid-crystalline polymer films 
with low energetic disorder. 

When spatial disorder is taken into account, in the 
form of randomly varying inter-chain distances, the ob- 
served crossover from increasing to decreasing mobility 
with temperature is no longer present. Instead, the 
mobility always increases with T, as was seen in the 
experimentscj. Because of the presence of a distribution 
of the relative position vectors, the effect of temperature 
is now enhanced, even though the mobility is still only 
weakly dependent on the electric field . Similar behaviour 
was also observed in films in which the chains deviate 
slightly from perfect nematic alignment perpendicular to 
the direction of the electric field. 

Finally, our numerical simulations in liquid-crystalline 



polymer films have helped us to distinguish between the 
effects of energetic and spatial disorder in such systems. 
Low energetic disorder, thus, high chemical purity of the 
chains, is responsible for the weak dependence of the 
mobility in the field. The inclusion of spatial disorder, 
associated with irregularities in the arrangement of the 
polymer chains within the film, does not make this de- 
pendence any stronger so long as the energetic disorder 
remains low. The influence of spatial disorder appears 
at the quantitative behaviour of the mobility, as it re- 
sults in decreasing considerably the absolute value of ln/x 



for a given range of fields and at a given temperature. 
We should note, however, that our model does not al- 
low the calculation of the exact value of the mobility, or 
make predictions about how many orders of magnitude 
the mobility might increase or decrease. 
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FIG. 1: Current transient for energetic disorder a = 2, with E = 3 x 10 5 V/cm, and T = 300K. The transit time tr is indicated 
by the arrow. 
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FIG. 2: Current transient for T — 300K and E = 3 x 10 5 V/cm. Curve (a) corresponds to the case of no energetic disorder 
(5e — 0), and curve (b) corresponds to a — 2. Curve (a) has been shifted along the y axis for clarity. 
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FIG. 3: The logarithm of the mobility ln(/i) as a function of yfE and for various temperatures T. The energetic disorder 
a = 1. /i is in arbitrary units. 
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FIG. 4: The logarithm of the mobility ln(^t) as a function of VE and for various temperatures T. The energetic disorder 
a = 2. /i is in arbitrary units. 
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FIG. 5: The logarithm of the mobility ln(/i) as a function of T. The circles correspond to a — 1, and the squares to a — 2. 
The filled symbols are for E — 2.6 x 10 V/cm, and the open symbols are for E — 4.6 x 10 J V/cm. fi is in arbitrary units. 
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FIG. 6: The logarithm of the mobility fi as a function of yfE and for various temperatures T for irregularly spaced polymer 
chains. The lines are a guide to the eye. The energetic disorder is equal to a = 2. //is in arbitrary units. 
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FIG. 7: The logarithm of the mobility fi as a function of \AE for irregularly spaced polymer chains, for T = 300K. The diamonds 
correspond to a = 2, and the circles to a = 3. The lines are a guide to the eye. ft is in arbitrary units. 
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FIG. 8: The logarithm of the mobility fi as a function of Vc and for various temperatures T, for randomly oriented and 
irregularly spaced polymer chains. The energetic disorder is equal to a = 2. The lines are a guide to the eye. ft is in arbitrary 
units. 



